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1. Introduction 

The quantisation of the ten-dimensional superstring using pure spinors as world-sheet 
ghosts [|lj has overcome many difficulties encountered in the Green-Schwarz (GS) and 
Ramond-Neveu-Schwarz (RNS) formalisms. Most notably, by maintaining manifest space- 
time supersymmetry, the pure spinor formalism has yielded super- Poincare covariant multi- 
loop amplitudes, leading to new insights into perturbative finiteness of superstring theory 




Counting fermionic zero modes is a powerful technique in the computation of loop 
amplitudes in the pure spinor formalism and has for example been used to show that 
at least four external states are needed for a non- vanishing massless loop amplitude Q. 
Furthermore, the structure of massless four-point amplitudes is relatively simple because all 
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fermionic worldsheet variables contribute only through their zero modes. In the expressions 
derived for the one-loop || and two-loop § amplitudes, supersymmetry was kept manifest 
by expressing the kinematic factors as integrals over pure spinor superspace Q involving 
three pure spinors A and five fermionic superspace coordinates 9, 

ffi-ioap = ((AA)(A7 m W0(A 7 n W0^ m „} , (1 x) 

^2-ioop = <(A 7 mnwr A)(A 7 s ^)^ m „ T pq T rs ) , 

where the pure spinor superspace integration is denoted by (. . . ), and A a (x,9), W a (x,6) 
and J r mn (x,9) are the superfields of ten-dimensional Yang- Mills theory. 



The kinematic factors in (1.1) have been explicitly evaluated for Neveu-Schwarz states 
at two loops P] and one loop 0, and were found to match the amplitudes derived in 
the RNS formalism ||]. This provided important consistency checks in establishing the 
validity of the pure spinor amplitude prescriptions. (Related one-loop calculations had 
been reported in |J.) 

In this paper, it will be shown how to compute the kinematic factors in ( |1.1| ) when 
the superfields are allowed to contribute fermionic fields, as is relevant for the scattering 
of fermionic closed string states as well as Ramond/Ramond bosons. It turns out that 
the calculation of fermionic amplitudes presents no additional difficulties, making ( |1 . 1| ) 
a good practical starting point for the computation of four-point loop amplitudes in a 
unified fashion. This practical aspect of the supersymmetric pure spinor amplitudes was 
also emphasised in |l(J, where the tree-level amplitudes were used to construct the fermion 
and Ramond/Ramond form contributions to the four-point effective action of the type II 
theories. 

This paper is organised as follows. In section 2, different methods to compute pure 
spinor superspace integrals are explored. These methods are then applied to the explicit 
evaluation of the kinematic factors of massless four-point amplitudes at the one-loop level 
in section 3, and at the two- loop level in section 4. In both these sections, the bosonic calcu- 
lations are briefly reviewed before separately considering the cases of two and four Ramond 
states. Particular attention will be paid to the constraints imposed by simple exchange 
symmetries. An appendix contains algorithms which were used to reduce intermediate 
expressions encountered in the amplitude calculations to a canonical form. 

2. Zero mode integration 

The calculation of scattering amplitudes in the pure spinor formalism leads to integrals over 
zero modes of the fermionic worldsheet variables A and 9. Both 9 and A are 16-component 
Weyl spinors, the A are commuting and the 9 anticommuting, and A is subject to the pure 
spinor constraint (A7 m A) = 0. The amplitude prescriptions |l], require three zero modes 
of A and five zero modes of 9 to be present, and a Lorentz covariant object 

j-a/3'y,5i...S B = (\ a \P\7Q Sl 9 Ss ) = T^ "^'^ 1 '" S ^ (2 1) 

was constructed such that the Yang-Mills antighost vertex operator 

V = (\7 m 9)(\ 7 n 9)(\Y9)(9 lmnp 9) has (V) = 1 . (2.2) 
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In this section, different methods of computing such "pure superspace integrals" are ex- 
plored. As an example, a typical correlator encountered in the two-loop calculations of 
section 4 is considered: 

F(ki, Ui ) = A;>^A;3A: r 4 ((A7 mnp9[r 'A)(A7 s ^i)(e7n a ^)(^2)(^3)(^7sn4)> (2.3) 

Here, k l and Ui are the momenta and spinor wavefunctions of the four external particles. 

2.1 Symmetry considerations and tensorial formulae 

One systematic approach to evaluate the zero mode integrals is to find expressions for all 
tensors that can be formed from (|2.1|). By Fierz transformations, one can always write the 
product of two 9 spinors as (9j^9), where 7^ denotes the antisymmetrised product of k 
gamma matrices. Due to the pure spinor constraint, the only bilinear in A is ^7^ A), and 
it is thus sufficient to consider the three cases 

((a 7 [5] a)(a{ 7 [1] or 7 ra or ^}9)(9j l3] 9)(9j l3] 9)) . (2.4) 

Lorentz invariance then implies that it must be possible to express these tensors as sums of 
suitably symmetrised products of metric tensors, resulting in a parity-even expression, plus 
a parity-odd part made up from terms which in addition contain an epsilon tensor. The 
parity-even parts may be constructed || starting from the most general ansatz compatible 
with the symmetries of the correlator and then using spinor identities along with the 
normalisation ( |2.2[ ) to determine all coefficients in the ansatz. Duality properties of the 
spinor bilinears can be used to determine the parity-odd part An extensive (and almost 
exhaustive) list of correlators is found in [11|, including the (Xj^9) and (Xj^9) cases of 
the above list: 

( (\'j mnpqr \)(\j u 9) {e^ 9h e){e^ kl e) ) 



1 / jrmnpqr . J_ mnpqr _ \ :mn trpq { srr srh ■ rr ri rr srh^ 

35 \°mnpqr 5! e mnpqr) 0f g Oj k {OiO u -f 0^0 U U I 



(2.5) 

[fgh]\jkl] 



' {\j mnpqr \)(\~f stu e) {9j f9h 6) (9j jkl 9) ) 



24 (jrmnpqr , J_ mnpqr \ cm srnp xqxt/srrxk _ srkrr^. 

' 35 \ u mnpqf ' 5! fc mnpqr) "j u f g u s u l \ u h u u u h u u> 



(2.6) 

[fgh][jkl](fgh^jkl) 



(Here, the brackets (fgh *-* jkl) denote symmetrisation under simultaneous interchange of 
fgh with ijk, with weight one.) The remaining correlator with the (Xj^9) factor can be 
derived in the same way, using an ansatz consisting of six parity-even structures. Taking 
a trace between the two 7^ factors and noting that 

Var{{\l mnpqr \){\l abcde 9) ...) = -4((X-f mnp ^ b X){Xj cde ^9) . . . > , 

one finds a relation to (|2.6|). This is sufficient to determine all coefficients in the ansatz, 
and the result is 

^ Xl mn ^ r X)(X^9)(9^ h 9)(9^ kl 9)) = f (C^>- + ^e mnpqr mnpqf) 

€£ PS i 6d 9 S k(.- S hti + 2 TO + C^fa&h^i - 3«)1 ( 2 - 7 ) 

. 3 9 k 19 \[abcde][fgh}[jkl\{fgh^jkl) 
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One may find it surprising that the derivation of these tensorial expressions only made 
use of properties of (pure) spinors, and of the normalisation condition fl2.2| ). However, it 
can be seen from representation theory that the correlator ( [2.1] ) is uniquely characterised, 
up to normalisation, by its symmetry. To see this, note that the spinor products A 3 
and 9 5 transform in 

A (a A /3 A 7) . Sym 3 s + = [ 0000 3] [ 100 oi] 

e^ 1 . . . e Ss] : Ait 5 s + = [00030] e [11010] . ^ 

(Here, A and 8 are taken to be in the S + irrep of SO(l,9), with Dynkin label [00001].) The 
tensor product of these contains only one copy of the trivial representation. This applies 
to any spinors A, which means that the pure spinor property cannot be essential to the 
derivation of the tensorial identities. The use of the pure spinor constraint merely allows 
for simpler derivations of the same identities. 

As an illustration of this approach, consider the correlator of eq. ( |2.3[ ). Leaving the 
momenta aside for the moment by setting F = k^k^kpkfF, the task is to compute 

p = ((A7 mnM[r A)(A7^i)(07n a ^)^7^2)(^3)(^^4)> . 
After applying two Fierz transformations, 

F = ( ^ <(A7™ [r| A)(A7^)(^7n a ^)(^ M ^)> (ui7 W 7c7^ 2 ) 

+ 3T16 ((^ mnm[rl mi Cde 0)(07n ab 8)(9^ kl e)) (u^JcdeW) 

+ OT <(A7 mnp9[r| A)(A7 cde/ ^)^7n a ^)(^ H ^)> (u^hcdefgW] 



one obtains a combination of the fundamental correlators listed in (^^), ( |2.6| ) and (p.7|). 
A reliable evaluation of the numerous index symmetrisations is made possible by the use of 
a computer algebra program. In doing these calculations with Mathematica, an essential 
tool is the GAMMA package |13[] , expanding the products of gamma matrices in a basis. 
The result consists of two parts, F = F^ + F^ e \ with 

F(6) = 5M( w i7 mpr ^)(u 3 7V) + 7io^C(«l7^2)(«37i^) + • • • 

-im^ aill2u ^ u ^ mpr ^A) (92 terms) (2.9) 

F{£) = -T2ok)O e n... J7 mpr (^i7 n --- i ^2)(n37V) + ... 

" W^ mpr ii... 4 6(^i7 i3 --- l9 ^2)(^37 nt2 ^ 8 i 9 n4) (34 terms) (2.10) 

The epsilon tensors in the second part can be eliminated using the fact that the Ui are 
chiral spinors: If all the indices on j^Ui are contracted into an epsilon tensor, one uses 

1 , 

. jk„,J"i _ ( \ 

'i\...iyl\---3k 

where 7 11 = ■^£i ...i 9 'y 10 '" 19 ■ More generally, if all but r indices of y^Ui are contracted, 



fc?.i...i,_/ii...i % 7 7 K - ll\...ly ■> \ Z - LL ) 



. . . . ^pi...vrh-ik~n (_\\KM)u ^L! ap—pw ■/, (2 12) 
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The result of these manipulations is 



P(£) = ~ 5Bo(«l7 mP ^2)(^37^4) " r±5 P r ( Ul ^ mi U 2 )(u 3li U A ) + 
11200* 



+ iim( u ^ ill2l3u 2)(^l ampr il i 2 i 3 u i ) (53 terms) (2.13) 



(Note that while the epsilon terms in the basic correlator formulae were easily obtained from 
the delta terms by using Poincare duality, this cannot be done here in any obvious way.) 
The last step in the evaluation of ( |2.3j ) is to contract with the momenta, F = k^k^kpkfF, 
and to simplify the expressions using the on-shell identities £V fc, = 0, kf = 0, = 0. It 
is shown in appendix [A. 2 that there are only ten independent scalars, denoted by B\ . . . B\o, 
that can be formed from four momenta and the four spinors u\ . . .114. With respect to this 
basis, the result is 

F = 48-10080 ( 695s i2(«i^3«2)(w3^i«4) H h 233si 3 {ui^ a u 2 ) («37a«4)) (7 terms) 

= Wmm ( 695 ' 775 > °> " 80 ' 356 ' 356 ' °' 233 ' 233 ' °)Bi -Bio - 
F(£) = 48Tot(- 23 '- 7 ' '- 16 ' 28 ' 28 '°' 7 ' 7 '°)Bi...Sio, 

F = ToS8o( 14 ' 16 ' '- 2 ' 8 ' 8 ' ' 5 ' 5 '°)Bi-S 1 o ) (2-14) 

where Sij = ki ■ kj. 

2.2 A spinorial formula 

While the derivation of tensorial identities for correlators of the form fl2.4|) is relatively 
straightforward and elegant, it may be a tedious task to transform the expressions encoun- 
tered in amplitude calculations to match this pattern. As seen in the example calculated 
above, this is particularly true if additional spinors are involved, making it necessary to ap- 
ply Fierz transformations. It is therefore desirable to use a covariant correlator expression 
with open spinor indices. Such an expression was given in M, 

fapyA-SB =N~ 1 \^ m ) aSl ^ n ) m (i P y S3 {j mn p) Si6B ] , (2.15) 

L J (a/3"/)[S 1 ...S 5 ] 

where N is a normalisation constant and the brackets ()[] denote (anti-)symmetrisation 
with weight one. (Note that the right hand side is automatically gamma-matrix traceless: 
any gamma-trace 

(Y)af3 X ( 7 m ) a[5ll (7 n ) /3|52l (7 P ) 7|53 (7mn P ) M5] = " ) ^ {imnp)^ (Y)^ 1 = 

vanishes due to the double-trace identity (7 a 6#) a (#7 af>c #) = 0, which follows from the 
fact that the tensor product (Alt 3 S + ) (g) S~ does not contain a vector representation and 
therefore the vector (V ; 7afe^)(^7 af>c ^) has to vanish for all spinors ip, and can also be shown 
by applying a Fierz transformation.) This prescription was originally motivated ^ by the 
fermionic expansion of the Yang-Mills antighost vertex operator V, 

V = T a/37 , 5l ... (55 A a A / W 1 ...0 s * (2.16) 

with 

T a fa,5i...8 6 = [(l m )aS 1 h n )pS 2 (Y)-/S 3 hrnnp)8 4 8 5 ] ia/3l)[Sl _ S5] , 
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where T is related to T by a parity transformation, up to the overall constant N. (Since 
T is uniquely determined by its symmetries, any covariant expression will be proportional 
to T, after symmetrisation of the spinor indices, and this is merely the simplest choice.) 

Equation ( |2.15| ) immediately yields an algorithm to convert any correlator into traces 
of gamma matrices or, if additional spinors are involved, bilinears in those spinors. It 
is, however, already very tiresome to determine the normalisation constant N by hand. 
The main advantage of this approach is that it clearly lends itself to implementation on 
a computer algebra system, which can easily carry out the spinor index symmetrisations, 
simplify the gamma products (again using the GAMMA package), and compute the traces. 
For example, 

L J (a/57) [01... 05] 

= ~h ^(7,7 m ) Tr( 7y7 ") Trfrrf) Tr( 7 ^ 7pnm ) + . . . 

- i TY( 7z7pnm7 ^ 7 - 7x7 m 73/ f>) (60 terms) 
= 5160960 . 

The correct normalisation is therefore obtained by setting N = 5160960. 

Returning to the example correlator ( p.3[ ), one finds that the calculation is by far 
simpler than with the previous method. After carrying out the symmetrisations (aj3j)[8i], 
one obtains 

NF=± Tr( 7x7 ab n7 , 7 ™[H )(u37g7 ^ 7sU4)(ni7 | S ] 7z7;)U2) + , . . 

-^(n 27b7 ^ 7g n 3 )(n l7s7?/7 ab n7:E7 ™[ r 727 s ] U4 ) , (24 terms) 

where elementary index re-sorting has reduced the number of terms from 60 to 24. Ex- 
panding the gamma products leads to 

NF = ^5 p r ( Ull m u 4 )(u 2l a U3 ) + ■■■ + ^(n l7 a4li2i3i % 2 )(n3 7 mpr ili2i3i 4«4) , (294 terms) 

which, in contrast to ( 12.10 ), contains no epsilon terms as there are not enough free indices 
present. Note that this intermediate result contains terms with with u\ paired with 113 
or U4, so it is not possible to directly compare to eqs. ( |2.9D and ( |2.13|) . However, after 
contracting with the momenta k^k^k^kf and decomposing the result in the basis B\ . . . Biq, 
one again obtains 

F = 10S8O ( 14 ' 16 ' °» - 2 . 8 ' 8 ' °» 5 ' 5 ' °)Bi ...Bio , (2-17) 
in agreement with ( p. 14 ). 

The algorithm just outlined will be the method of choice for all correlator calculations 
in the later sections of this paper and can easily be applied to a wider range of problems. 
The only limitation is that the larger the number of gamma matrices and open indices 
of the correlator, the slower the computer evaluation will be. For example, the correlator 
considered in eq. (5.2) of JTT|], 



m4?i4 , 



t mnmim . . .m 4 n 4 _ ( (y^im ^ {Xf*f*"*0) (Af f" 3 ™ 3 6) (<9 7 m 7 n 7p(?r7 r ' 

2 I vanA.m\n\...m^n^ 1 m,nm\n\...va^n^\ in -i o\ 

— — 45 VI l 8 2 £ ) ' \ Z - L °) 

can still be verified with this method but this already requires substantial runtime. 
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2.3 Component-based approach 

A third method to evaluate the zero mode integrals consists of choosing a gamma matrix 
representation, expanding the integrand as a polynomial in spinor components, and then 
applying ( 2.15| ) to the individual monomials. This procedure seems particularly appealing 



if at some stage of the calculation one works with a matrix representation anyhow, in 
order to reduce the results to a canonical form (e.g. as outlined in appendix [A|). An 
efficient decomposition algorithm (of k u\U2U^Ui scalars, say) only needs a few non-zero 
momentum and spinor wavefunction components to distinguish all independent scalars, and 
therefore k and u can be replaced by sparse vectors. Furthermore, a trivial observation 
allows for a much quicker numeric evaluation of correlator components than a naive use 
of (gig) : In view of (|I|) , one can equivalently compute the components of the parity- 
transformed expression V = (Xj m 9)(X^ n 9)(X^f p 9)(6'y rnn p8), where A and 9 are spinors of 
chirality opposite to that of A, 9. In the representation given in appendix [BL V coincides 
with V\ x ^ g _ s , and 

V = 192 \ 9 \ 9 \ 9 9 1 9 2 9 3 9 A 9 9 + • • • + 480 A^A 3 ^ W 3 # 15 + . . . (100352 terms) 

The monomials in the fermionic expansion of V then correspond to the arguments of 
non-zero correlators, and the coefficients of those monomials are, up to normalisation and 
symmetry factors, the correlator values. 

Unfortunately, it turns out that the complexity of typical correlators (e.g. the one 
given in ([2.3;)) makes it difficult to carry out the expansion in fermionic components in 
any straightforward way and limits this method to special applications. For example, the 



coefficients in ( 2.18 ) can be checked relatively easily by choosing particular index values, 
such as 

<(A7 P 7 1 ^)(A7 g 7 21 ^)(A7 r 7 34 ^)(^7 7 7 W r7 43 ^)> 

= (12 X 1 X 1 X 1 6 1 6 9 8 w 9 n 9 12 + ■■■ + 12 A 16 A 16 A 16 0W0 8 16 ) = ^ . 

(For fixed values of pgr, one gets no more than about 10 5 monomials of the form A 3 # 5 ). 
This approach may thus still be helpful in situations where the result has been narrowed 
down to a simple ansatz. 

3. One-loop amplitudes 

The amplitude for the scattering of four massless states of the type IIB superstring was 
computed in the pure spinor formalism as 

A = KK J J £ Z2 I d i Z3 I fPztYlGfazj)*-**, (3.1) 

where G(zi,Zj) is the scalar Green's function, and the kinematic factor is given by the 
product KK of left- and right-moving open superstring expressions, 

ifi-ioop = {(\A l )(\ 1 m W 2 )(\ 1 n W^ m n) + (cycl(234)) . (3.2) 
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Here the indices 1 ... 4 label the external states and "• • ■ + (cycl(234))" denotes the addition 
of two other terms obtained by cyclic permutation of the indices 234. The spinor super- 
field A a and its supercovariant derivatives, the vector gauge superfield A m = \"im ] D a Ap 
as well as the spinor and vector field strengths W a = jQ(y m ) al3 (DpA m — d m Ap) and 
Fmn = \{"imn) a pDaW^ = 2<9[ m j4 n ] , describe ten-dimensional super- Yang-Mills theory. 
The physical fields of this theory, a gauge boson and a gaugino, are found in the leading 
components A m \ = £ m and W a \ = u a and correspond to the Neveu-Schwarz and Ramond 
superstring states. 

The superfields A a and W a as well as the gaugino field u a are anticommuting. 1 To 
facilitate computer calculations involving polynomials in the spinor components, and for 
easier comparison with the literature, it will be more convenient to work with commuting 
fermion wavefunctions u a . Fortunately, as the kinematic factors with fermionic external 
states are multilinear functions of the distinctly labelled spinors Uj, it is straightforward to 
translate between the two conventions: Any monomial expression in u\ . . . U4 (and possibly 
fermionic coordinates 0) corresponds to the same expression in u\ . . . U4, multiplied by the 
signature of the permutation sorting the iii (and any 9 variables) into some fixed order, 
such as (6---e)u^u^ 2 u^u2 4 . 

Choosing a gauge where a A a = 0, the on-shell identities 

2D (a A p) = 7 ™jA m , D a W P = \{l mn )jFmn 
have been used to derive recursive relations jR], [l4], 15 1 for the fermionic expansion 



4 W) = ^ m 0)aA^ , A (n) = I ( 7wW Kn-l)) w a(n) = .l^n^^^-!) 

where = ^6 an • • • 6 ai (D ai ■ ■ ■ D an f)\. These recursion relations were explicitly solved 
in |]10| , reducing the fermionic expansion to a simple repeated application of the derivative 
operator m « = \{e lm ^9)d p : 

A {2k) _ 1 unk] q/- 
■^rn — (2k)\ IV m , 



With this solution at hand, one has all ingredients to evaluate the kinematic factor ( p.2| ) 
for the three cases of zero, two, or four fermionic states. 

3.1 Review: four bosons 

The kinematic factor involving four bosons was considered in and this calculation will 
now be reviewed briefly. First, note that the outcome is not fixed by symmetry: The result 
must be gauge invariant || and therefore expressible in terms of the field strengths F\ . . .F^. 
The cyclic symmetrisation in ( |3.2|) yields expressions symmetric in F2, F3, F4, and acting 
on scalars constructed from the Fi only, the (234) symmetrisation is equivalent to complete 
symmetrisation in all labels (1234). Thus the result must be a linear combination of the 



1 Thanks to Carlos Mafra for pointing this out. 
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two gauge invariant symmetric F 4 scalars, namely the single trace Tr (i^n F2F3F4)) and 
double trace Tr (.Fni*^) Tr(F 3 F^\), leaving one relative coefficient to be determined. 

Since all four states are of the same kind, one may first evaluate the correlator for one 
labelling and then carry out the cyclic symmetrisation: 

Kifol = ((^i)(\j m W 2 )(\ 7 n W 3 )F 4 , mn )\ iB + (cycl (234)) . 
The different ways to saturate 9 5 result in a sum of terms of the form 

Xabcd = ((XA^)(Xj m W^)(\j n wP)^X) (3-4) 
with A + B + C + D = 5 and A, B, C odd, D even: 

<(A^i)(A7 m W 2 )(A7 n ^ 3 )^4, m n>| 4B = X 3110 + X 13W + X U30 + X 1U2 . 

Note that X\ 3 iq and X1130 are related by exchange of the labels 2 and 3. This exchange 
can be carried out after computing the correlator, an operation which will in the following 
be denoted by 7r 23 . Using (3.3) for the superfield expansions and replacing d m — > ik m , one 
obtains 

X3110 = -isFi n F^ s Ft u X 3110 , X3110 = ((A 7 [S p ^)(A7 |ul 7 r ^)(A7a^(^ am ^)> , 
Xni2 = -^ikiCFp g F?sFiuXiU2 , X im = <(A 7 [m| 7 P ^)(A7 |al 7 r ^)(A 7 ^)(^^)> , 
*i3io = -^kUlF^XuXi^ , X 13W = <(A 7 [t| 7 m ^)(A7 |u] 7 r ^)(A7^)( • 



The method outlined in section 2.2 is readily applicable to these correlators. For example, 
for X 3 m, the trace evaluation yields 

^3110 



X 3 mo — N 1 



± Tr( 7 a7 2 ) Tr( lxyzl a ™) TY^^ 1 ) Tr( 7 ^ r7 M ) + ■ ■ • 
• • • + i Tr( 7 ["l 7rs72y:EW l i ] 7 a; 7a7 i/ 7 m " a 7 z )] (60 terms) 

— ( 2^K m PS n 1 l_X mn XPq ]_xmnt-pq , 26 rmnri)s\ 

\ 3 5°rs °tu 3 15 u tu u rs 45°rs °tw 315 V °tu) [mn]\pq][ra][tu](pg^rs) 

Upon contracting with the field strengths, momenta and polarisations, and symmetrising 
over the cyclic permutations (234) (with weight 3), one finds that all three contributions 
are separately gauge invariant: 

X 3 no + (cycl(234)) = Tr(F (1 F 2 F 3 F 4) ) + ^ TV(F (1 F 2 ) TY(F 3 F 4) ) 

Xnu + (cycl(234)) = Tr(F (1 F 2 F 3 F 4) ) + ^ Tr(F (1 F 2 ) Tr(F 3 F 4) ) 

(1 + vr 2 3)X 13 io + (cycl(234)) = (4Tr(F (1 F 2 F 3 F 4) ) - Tr(F (1 F 2 ) Tr(F 3 F 4) )) 

The sum X3110 + -X"iii2 has the right ratio of single- and double-trace terms to be propor- 
tional to the well-known result t^F 4 , and the last line exhibits the right ratio by itself. The 
overall kinematic factor is therefore 

^i-ioop = "2^0 (4TV(F (1 F 2 F 3 F 4) ) -Tr(F (1 F 2 )Tr(F 3 F 4) )) = -j^t s F\ (3.5) 

in agreement with the expressions derived in the RNS and Green-Schwarz [ jl7f | for- 
malisms. 
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3.2 Four fermions 

The four-fermion kinematic factor could be evaluated in the same way as in the four-boson 
case by summing up all terms Xabcd, A + B + C + D = 5, now with A, B, C even 
and D odd. Note however that this time, the outcome is fixed by symmetry: The cyclic 
symmetrisation in fl3.2j ) leads to a completely symmetric dependence on u 2 , U3, U4, and 
therefore to a completely antisymmetric dependence on 112, 113, 114. Acting on scalars of 
the form k 2 u\U2Uj,u^^ antisymmetrising over [234] is equivalent to antisymmetrising over 
[1234], and there is only one completely antisymmetric k 2 uiU2Usu^ scalar. Without further 
calculation, one can infer that the kinematic factor is proportional to that scalar, 

which of course agrees with the RNS amplitude (see e.g. fl(|, eq. (3.67)). 

3.3 Two bosons, two fermions 



In evaluating ( |3.2| ) for two bosons and two fermions, the cyclic symmetrisations affect 
whether the W and T superfields contribute bosons or fermions. Only the label of the A a 
superfield stays unaffected, and one has to choose whether it should contribute a boson 
or a fermion. Since its fermionic expansion starts with the bosonic polarisation vector, 
A\. a ~ (^i^)aj the calculation can be simplified by choosing a labelling where particle 1 is 
a fermion. (Of course, the final result must be independent of this choice.) The assignment 
of the other three labels is then irrelevant and will be chosen as /i/2^3^4- Writing out 
the cyclic permutations, two of the three terms are essentially the same because they are 
related by interchange of the labels 3 and 4. The kinematic factor is then 

^K(/l/2« = (1 + 7r 34 )<(A4 eVen) )(A 7 ^rt 

+ ((A^S even) )(A 7 m ^3 (odd) )(A7 n ^i odd V£n ) > • 

Unlike in the four-fermion calculation, the result is not fixed by symmetry. There are five 
independent kui^F^Fi scalars (see appendix [A], eq. ( A.6[ )), denoted by C\ . . . C5, and there 



are two independent combinations of these scalars with the required [12] (34) symmetry. 
Expanding the superfields and collecting terms with 6 5 , the first line yields a combination 
of terms Xabcd with A, B, D odd and C even. There is only one 9 5 combination coming 
from the second line, which will be denoted by X' 2111 = (— ^24)^2111 : 

-^Woop = (1 + 7r 34) (XlOlO + -^2210 + -^2030 + -^2012) + X' 2 i U , 



with the correlators 



-^4010 


_ i ul r,3 z-3 h.4 z-4 v 

— gQt g ^C, c 'v ?n l, n ^-4010 > 


-^4010 


-^2210 


= ~ T2 k b k dCe k mCn- X 2210 , 


-^2210 


AT2030 


= ~iE k b k dCe k mCt X 2030 , 


-^2030 


-^2012 


= ~ T2 k b(c k m k dCt X 2012 , 


-^2012 


^2111 


— 6 b^c^d^e h 'm^2\\\ 1 


^2111 



((X^ a 9)(9 la p W)(e 7pUl )(X^ m u2)(Xj n h bc 0)) 
<(A7^)(^i)(A7 [m| 7^)(^7 c n2)(A 7 |n] 7 d ^)) 

< (Xj a e) (e laUl )(\^ m u2)(\7 n] j bc o)(Ojc de e)) 

<(A 7 ^)(^i)(A7 [m n 2 )(A 7 nl 7^)(07n d ^)> 

((\7 a e)(e laUl )(x^ m h bc e)(x^ n h de o)(e ln u2)[ 
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(The numerical coefficient in AT 2111 includes a sign coming from the 0, u ordering: there 
is an odd number of 9s between u\ and u 2 .) Evaluating these expressions as outlined in 
section ^2, the spinor wavefunctions U{ present no complication. The last part takes the 
simplest form: One finds 

((A 7 ^)(0 7a n 1 )(A 7 m 7 ^)(A7 n 7 d ^)(^ 2 )> = - ^ (2<t[> l7e p 2 ) + t( Ull c ^u 2 )) 
and therefore 

XL 



L 2111 

The result for X4010 is 



480 



l ^( Ul7 VV>+^l7VV 



X.. 



4010 



k 5 mn(uil C U2) ~ m$mq(ui7 n U 2 ) + 4n 5% n (Ul7% 2 ) - aggoCC"!^^) 



360 1 
1 xb 



bcmnq 



U2, 



[be] [mil] 

For the evaluation of X221O) it is useful to consider the more general correlator 



'(A7^)(0 7 ani)(A 7 [ m l 7 ^)(A7 |n] 7 d ^)(^7^2) 



+ 



11 rm / 



bedemnx 



201600 °z V"17 2/ ' 403200V" 1 ' 

9676800 ^bcdemniii2i3i4 (^17 1 2 3 4 U 2 ) 



«2)] 



(27 terms) 



+ 



[mra] [fee] [de] 

^-bedemnxixi?. 



2419200 



This time, even using the method of section 2.2, there are sufficiently many open indices 
and long enough traces for epsilon tensors to appear. Using eqs. (2.11) and ( |2.12| ), they 
can be re-written into 7 [ 5 ' 7 ] terms: 



;(A 7 ^)(^i)(A7 [m| 7^)(A 7 |nl 7' i ^)^7^2) 



+ 



1 xm 



5™( Ull bcden u 2 ) + 



16800 "x 



33600 



(uii 



bedemnx 



U 2 , 



[inn] [be] [de] 



(27 terms) 



A good check on the sign of the epsilon contributions is that A^iii is recovered when 
contracting with rj nx , involving a cancellation of all 7® terms. To obtain -X221O5 one 
multiplies by — r/ cx : 



X- 



2210 



720^(^7^2) + mo^ n ( Ul j e u 2 ) + ^f^T"^) + 20I60 4 (^17 6 "'" "2 i 



1 1 rfe / den \ 1 11 At emn \ . 1 i 
+ 2880 d m(«l7 U 2 ) + 2 Qm d d( u H ^2) + l48o(Ul7 



bdemn„ 



u 2. 



[de] [mri] 

For the calculation of AT2030 an d -^2012 j one may first evaluate a more general correlator 
({\^ a Q)(Q^aUi)(\^ m u 2 )(\i a ^ bc 6)(6"f x ^ de 6)) and then contract with r] cx and rj nx , respec- 
tively. The results are 



X 



2030 



-72oCn(^i7%) + 3sA(«l7 e «2) " TS0^K7 n n 2 ) - ToTM^K7 &e ^2) 
■^^(tiiT^tia) - -^(u l7 emn « 2 ) + ^o(n l7 6de "" 



u 2, 



[ran] [de] 



x. 



2012 



2M^m( n l7^2) + 2m5L(«17^2) " THo^^ 7 "^) + 20I6 S i (^17^2 



11 
10080' 



<t(ui7 cc %) + jM^S b d ( Ull cem u 2 ) - ^-l Ull b ^ m u 2 ) 



[be] [de 
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After multiplication with the momenta and polarisations, all individual contributions are 



gauge invariant and can be expanded in the basis C\ . . . C5 listed in (A. 6): 



(1 + 7T34)X 4 oio = 4§^4o(-6, -16,-40,6,0) Cl ...c 5 
(1 + 7T34)X 2210 = 3g ^ 3g (-18,-104,-176,18,0) Cl ...c s 
(1 + 7r 34 )A 2 o3o = idko (-21, 42, -42, 21, Ota ...c 8 
(1 + vr 34 )A 2 oi2 = 48^4o(-39,78,-78,39,0) Cl ...C5 

^2111 = _ Tt|2o( 1 ' ' 4 ' ~ 1 ^)Ci...C 5 



The sum can be written as 

-2B2F 



^1-loop = X 2111 = _ 3Mo( 1 ' ' 4 ' - 1 )°)Ci...C 5 



1 

1920 



Sn(u 2 i 3 (^ 2 + h)ii U l) + s 23(«2i£ 4 (#2 + ^4)^3 U l) ) ( 3 - 6 ) 



and again agrees with the amplitude computed in the RNS result, see [16] eq. (3.37) 



4. Two-loop amplitudes 

The pure spinor formalism was used in |2| to compute the two-loop type-IIB amplitude 
involving four massless states, 

/ 4 r exp ( — X^j 7 ki ■ kj G(zi, Zj)) 
d 2 n n d 2 n 12 d 2 n 2 2l[j d 2 Zi ^ (detlmft) 5 - K 2-ioo P (ki,Zi) , 



where is the genus-two period matrix, and the integration over fermionic zero modes is 
encapsulated in 

^2-ioop = A 12 A 34 <(A7 mnMr A)(A 7 W 1 )^ 2 , mri ^3, TO ^4,rs> + (perm(1234)) (4.1) 
= Ai 2 A 34 ^i 2 + Ai 3 A 24 ^i 3 + Ai 4 A 23 i^i 4 . (4.2) 

The kinematic factors K12, if 13, i^i 4 are accompanied by the basic antisymmetric biholo- 
morphic 1-form A, which is related to a canonical basis u\,u 2 of holomorphic differentials 
via Aij = A(zi,Zj) = u>i(zi)u)2(zj) — uj2{z.i)uj\{zj). The superfields W" and Ti, mn are the 
spinor and vector field strengths of the i-th external state, as in section |3| One encounters 
superspace integrals of the form 

Y(abcd) = ((\7 mnpqr X)(X7 s W a )^ mn ^ pg ^ rs ) . (4.3) 

The symmetries of the A 3 combination jl in this correlator include the obvious symmetry 
under win *-* pq, and also {X^ mn Pi r X)(X-f s] )a = (this holds for pure spinor s A and can be 
seen by dualising, and holds for unconstrained spinors A as part of a A 3 # 5 scalar, as seen 
from the representation content (|2.8|)), and allow one to shuffle the T factors: 



Y(abcd) = Y(acbd) , Y(abcd) + Y(acdb) + Y(adbc) = . (4.4) 
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4.1 Review: four bosons 

The case of four Neveu-Schwarz states was considered in Q and will be briefly reviewed 
here. As all three kinematic factors If 12, K13 and K14 are equivalent, it is sufficient to 
consider K\ 2 in detail. With all external states being identical, the symmetrisations of 
( |4.1| ) can be carried out at the end of the calculation: 

Kg = 4 (W [1 f 21 f [3 f 41 )| 4B + 4{^ 4] f |1 f 2] }| 4B 

= (1 - 7Q 2 )(1 " vr 34 )(l + 7T137T24) (WlttsFi) | 4B 

Expanding the superfields and adopting the notation 

Y ABCD (abcd) = ((A 7 ™ r A)(A 7 ^( A ))^„4S^) , 

the Neveu-Schwarz states come from terms of the form Yabcd = ^a_bcd(1234) with A odd 
and B, C, D even. Using the shuffling identities fl4.4j) to simplify, one obtains 



= ^5000 + ^1400 + ^1040 + ^1004 + ^3200 + ^3020 + ^3002 + ^"l220 + ^1202 + ^1022 
= (1 + 7T 23 )(1 - 7r 24 ) + ^1400 + ^3200 + *1022) , 

and therefore Kf^ can be written as the image of a symmetrisation operator S^b'- 

^12 = ^4B (1^5000 + ^1400 + ^3200 + ^1022) 
5 4 B = (1 - VTi 2 )(l - vr 34 )(l + 7ri 3 7r 24 )(l + 7T 23 )(1 - 7T 24 ) 

It is worth noting at this point that, on the sixteen-dimensional space of Lorentz scalars 
built from the four field strengths and two momenta, the symmetriser £43 has rank four. 



The correlators were computed in J6|], using the method outlined in section 2A. Two are 
zero, I5000 = ^1400 = 0, and the remaining ones are 

I3200 = TkklF^klF^F* ((A 7 ™ r A)(A 7 s 7 a ^)(^ c ^)(^n e ^)>, 
^1022 = kFl h Fl n klF* d ktF* f ((A 7 ™[ r A)(A 7 s ] 7 a ^)(0 7g c ^)(^ 7s e ^)>. 

In reducing those two contributions to a set of independent scalars, one finds that they 
both are not just sums of (k ■ k)F 4 terms but also contain terms of the form k • F terms. 
The latter are projected out by the symmetriser S4B, and the result is 

Kff = S 4B (Y 3W0 + Y 1022 ) = ^(s l3 - s 23 ) (4Tr(F (1 F 2 F 3 F 4) ) - Tr(F (1 F 2 ) Tr(F 3 F 4) )) , 

= ^(SlS — S 23 )t8-F 4 • 

By trivial index exchange, one obtains K\% and ifu, and the total is 

K t B \oo? = m (( s i3 - s 23 )A 12 A 34 + (s 12 - s 23 )Ai 3 A 24 + (s 12 - s 13 )A 14 A 23 )t 8 F 4 , (4.5) 

a product of the completely symmetric one-loop kinematic factor tsF and a completely 
symmetric combination of the momenta and the A ^ . 
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4.2 Four fermions 

The calculation involving four Ramond states is very similar to the bosonic one. Focussing 
on the K\2 part, the symmetrisations in (^^) can again be rewritten as action of sym- 
metrisation operators on the correlator of superfields with one particular labelling: 



Kg(ui) = (1 - tt 12 )(1 - 7r 34 )(l + 7r 13 7r 24 ) (W^^i) 

= 4(1 - 7T 12 ) (Wi^ 3 ^4) 



The last step follows from the fact that all scalars of the form k^u 4 (see appendix A.2), 
and therefore all fc 4 « 4 scalars, are invariant under 7ri37r 24 and have 7Ti 2 = ^34. This time, 
on expanding the superfields, one collects the terms Yabcd with A even and B, C, D odd. 
After using ( [4.4[ ) to simplify, 

< W 1^8^4>| = Y21U + y ° 311 + y ° 131 + y ° 113 

= (1 + 7r 23 )(l - vr 24 ) (|l2in + Yani ) , 

and after translating to commuting wavefunctions ui, which multiplies every permutation 
operator with its signature, one obtains 

Kg{m) = S 4F (\Y 211x { Ui ) + ro3n(«i)) , S 4F = 4(1 + ^ 12 )(1 - 7r 23 )(l + ^ 24 ) . 

This symmetriser has rank three, and the result is again not determined by symmetry. 
Two correlators have to be computed: 

Y 2U i(u i ) = (-2)klk 2 m k 3 p kt ((A7 mnw[r A)(A7 sl 7 a ^)(^i)^7nU 2 )(07^ 3 )(^7^ 4 )> 
W(0 = (-l)k 2 a k 2 m k 3 p kt ((A7 mnw[r A)(A7 s ^i)(07n a ^)(^2)(^7^ 3 )(07^4)) 



With four fermions present, the method of section |2.2| is preferred as it does not involve re- 
arranging the fermions using Fierz identities. The first correlator was covered as an example 
in that section, and the second one can be evaluated in the same fashion. Expressed in the 
basis listed in ( A.5 ), the results are 

¥2111 = 5^0 (-19, -21, 21, 19,-17,-17,0, 0,0,0) Bl ... Bl0 , 
Y 3u(ui) = ^(-14,-16,0,2,-8,-8,0,-5,-5,0)^...^. 

After acting with the symmetriser 5 4 p, one obtains the same u 4 scalar encountered in the 
one-loop amplitude, 

Kg (m) = «S 4F ( 3^2111 (^) + Yo3ii(ui)) = ^(-1,-2,1, 2,-1, -2, 0,0,0,0) Bl ... Bl0 

= 55(323 - S13) ((«1^3 U 2)(^3^1^4) - (ui^2U 3 )(u2h v "i) + ("l h 11 "^ («2^1«3)) • 

The if 13 and K14 parts again follow by index exchange, and the total result 

K 2Aoop{ u i) = i(( s 23 - Si 3 )Ai 2 A 34 + (s 23 - Si 2 )Ai 3 A 24 + (s i3 - Si 2 )Ai 4 A 23 ) 

x {{u\^uz){u^ x U4) - (1*1^3) 0«2#iii4) + (uih u 4:){u 2 ^ 1 u 3 )) (4.6) 

is again a simple product of the one-loop kinematic factor and a combination of the Ajj 
and momenta. 
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4.3 Two bosons, two fermions 

As in the one-loop calculation of section |3.3| , in the mixed case one has to pay some attention 
to the permutations in Q4.1| ) since they affect which superfields contribute fermionic fields. 
The complete symmetrisation makes it irrelevant which labels are assigned to the two 
fermions, and the convention $\fobzb/± will be used here. The kinematic factor K 2 J^ 2F is then 
distinguished from the other two, Kff 2F and Kff 2F . Carrying out the symmetrisations in 
(|4.l| ) and using the identities (|4.4|) , one finds 



K 12 (U1,U2,(3,(4:) = (1 - tti 2 )(1 - 7r 34 )i^, 
ifl3(«l,«2,C3,C4) = (2 • 1 + 7Ti2 + vr 34 + 27ri27r 34 )K , 

Ku(ui,U2,(3,(4:) = (1 + 27T12 + 27T 3 4 + 7Ti2Vr 3 4)K , 

where, schematically, 

K = ( ( A 3 Ty i (CVen ' ) )^° dd ' ) jr( evcn ) j^r(even) ^ + ^^^(oddj^evenj^oddj^odd)^ _ 

In translating to commuting variables u\ and u 2 , the permutation operator ixyi changes 
sign, and therefore 2 

Ki 2 {ui,u 2 , Ca, U) = (1 + tti 2 )(1 - vr 34 )^ , 

^13(^1, "2,(3^4) = (2 • 1 - 7T12 + vr 34 - 27Ti 2 7r 3 4)if , 

Ku(lll,U2, C3, C4) = (1 - 27T12 + 27T 34 - TTuTTs^K . 

Expanding the superfields, the contributions to K are: 

^4ioo = -^klk l d k 2 m F^ s ((A7 mnp9r A)(A 7 s 7 a ^)(^7 c ^)(^7 c ni)(^7nn2)> 
^0500 = ^k 2 m k 2 a k 2 c F^ s ((A 7 ™ r A)(A 7 ^i)(^n a ^)(^ c ^)(^2)> 
^0140 = igk 2 m k 3 p k 3 a F? d F r \ <(A 7 ™ r A)(A7 s n 1 )(0 7n n2)(^ a ^)(^ c ^)> 
^0104 = <(A7 mw A)(A 7 ^ 1 )(0 7n n2)(^ 7 |/^)(0 7fe c ^)) 
^2300 = izklklklF^ ((A7 mnp ' r A)(A 7 s 7 a ^)(^i)(^^)(^«2)} 
I2120 = iklk^F^ ((A 7 m "^A)(A 7 s 7 a ^)(^n 1 )(0 7n n 2 )(^ c ^)> 
F2102 = <(A 7 mnwr A)(A 7 s 7 a ^)(0 7b n 1 )(^ 7n u 2 )(0 7 | s] c ^)) 

^0320 = UlklklF^ <(A 7 ™ f 'A)(A 7 s u 1 )(^ 7n a ^)(^2)(^ c ^)> 
W = iiklklF^F^ <(A 7 m "^A)(A 7 ^ 1 )(e 7n a ^)(^ 2 )(e 7 | s] c ^)) 
W = !^^^ 3 6 i^4| <(A 7 — A)(A 7 s « 1 )(0 7n u 2 )(0 7? °^)(0 7s] cd 0)> 
*3on = j-2 k l F cdFLklk 2 [r \ <(A 7 m "^A)(A 7 s 7 a ^)(^ c ^)(0 7c n 1 )(e 7n n 2 )> 
^1211 = §F a 3 b CF c 4 ^^f r| <(A 7 m "^A)(A 7 s 7 a ^)(0 7n c ^)(^n 1 )(^ 7 | sF2 )> 

n031 = i-2 F lAnk\k\k 2 A <(A 7 ™ r A)(A 7 S 7 a ^)(^ C ^)(^ 1 )(0 7 | sF 2)> 

^013 = j-2 F a b FLklk 2 c k 2 A ((A 7 m ^A)(A 7 s 7 a ^)(^u 1 )(0 7 | s] c ^)(0 7(i n2)> 

2 This sign change is crucial to avoid the erroneous conclusion that the two-boson, two-fermion kinematic 
factor cannot be of the same product form as in the four-boson or four-fermion cases, which would be in 
contradiction to the supersymmetric identities derived in mM. 
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These correlators can be evaluated exactly as described in section 3.3 . One finds that 
^0500 = ^0140 = ^0104 = 0, and the sum of the remaining terms reduces to 



K — I4100 + ^2300 + ^2120 + ^2102 + ^0320 + ^0302 + ^0122 + ^3011 + ^1211 + ^1031 + ^1013 

= 3gg(si2 + «is) x (1,0,4, -l,0) Cl ...c s • 
After applying the symmetrisation operators, 

(l + vr 12 )(l-7r 3 4)^ = i|o(si2 + 2si3) x (1, 0, 4, -1, 0) Cl ...c 5 , 
(2 • 1 - 7T12 + ^34 - 2^i27r 3 4)K = ^q(2s 12 + si 3 ) x (1,0,4, -l,0) Cl ...c 5 , 

(1 - 27T12 + 27T 3 4 - TT 12 TT U )K = j!^(s 12 ~ S 13 ) X (1, 0, 4, -1, 0) Cl ...C 5 , 

the total kinematic factor is seen to be 

K 2 -loop(ui,U 2 , C3, C4) = - ifo (( S 23 - S13) A12A34 + (s 2 3 ~ S12) Ai 2 A 34 + (s 13 - S 12 ) A12A34) 

x (1,0, 4,-1, 0) Cl ...c 5 (4.8) 
and displays the same simple product form as in the four-boson and four-fermion case. 

5. Discussion 

In this paper, different methods were discussed to efficiently evaluate the superspace inte- 
grals appearing in multiloop amplitudes derived in the pure spinor formalism. Extending 
previous calculations || restricted to Neveu-Schwarz states, it was then shown how the 
treatment of Ramond states poses no additional difficulties. 

While the bosonic calculations of |], [?j have, in conjunction with super symmetry, 
already established the equivalence of the massless four-point amplitudes derived in the 
pure spinor and RNS formalisms, it would be interesting to make contact between the 
results of sections PO / fO and two-loop amplitudes involving Ramond states as computed 



in the RNS formalism (see for example [19|). 

The assistance of a computer algebra system seems indispensible in explicitly evaluat- 
ing pure spinor superspace integrals. To avoid excessive use of custom-made algorithms, 
it would be desirable to implement these calculations in a wider computational framework 
particular adapted to field theory calculations |20[| . 

The methods outlined in this paper should be easily applicable to future higher-loop 
amplitude expressions derived from the pure spinor formalism, and, it is hoped, to other 
superspace integrals. 
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A. Reduction to kinematic bases 



In calculating scattering amplitudes one encounters kinematic factors which are Lorentz 
invariant polynomials in the momenta, polarisations and/or spinor wavefunctions of the 
scattered particles. It can be a non-trivial task to simplify such expressions, taking into 
account the on-shell identities ^ ki = 0, kf = 0, ki ■ d = 0, ^tij = 0, and, in the case of 
fermions, re-arrangements stemming from Fierz identities. 

More generally, one would like to know how many independent combinations of some 
given fields (subject to on-shell identitites) there are, and how to reduce an arbitrary expres- 
sion with respect to some chosen basis. This appendix outlines methods to address these 
problems, with an emphasis on algorithms which can easily be transferred to a computer 
algebra system. These methods are not limited to dealing with pure spinor calculations 
but the scope will be restricted to amplitudes of four massless vector or spinor particles in 
ten dimensions. 



A.l Four bosons 

It is not difficult to reduce polynomials in the momenta and polarisations to a canonical 
form. The momentum conservation constraint Yli = is solved by eliminating one 
momentum (for example foj), all kf are set to zero, and one of the two remaining quadratic 
combinations of momenta is eliminated (for example s 23 — > —s\ 2 — S13, where s^,- = ki ■ kj). 
Then all products ki ■ Q are set to zero, and one extra k ■ C product is replaced (when 
eliminating k±, the replacement is k 3 ■ £ 4 — > (— k\ — k 2 ) • £4)- The remaining monomials are 
then independent. (This is at least the case with the low powers of momenta encountered 
in the calculations of sections || and ||, where there are enough spatial directions for all 
momenta/polarisations to be linearly independent.) 

The implementation of these reduction rules on a computer is straightforward. The 
easiest way to obtain scalars which are also invariant under the gauge symmetry ki — * Q 
is to start with expressions constructed from the field strengths F^ b = 2d^ a $. For the 



one- loop calculations of section 3.1, the relevant basis consists of gauge invariant scalars 



containing only the four field strengths F\ . . . F4. One finds six independent combinations, 

Ti(F 1 F 2 F 3 F i ) Tr(F 1 F 2 )Tr(F 3 F 4 ) 
Tr(F x F 2 F 4 F 3 ) Tr{F 1 F 3 )Tr(F 2 F i ) 
Tr{F x F 3 F 2 F A ) TrGFxi^TrCFsFg) 



In the two-loop calculations of section ^lj, all monomials have two more momenta. There 
are sixteen independent gauge invariant scalars of the form kkFiF 2 F 3 F4, and twelve of 
them may be constructed from the previous basis by multiplication with s± 2 and S13: 
Ai = si 2 Tr(FiF 2 F 3 F4), A 2 = sx 3 Tr(FxF 2 F 3 F4,), etc. One choice for the additional four is 

A 13 = k 3 ■ F\ ■ F 2 ■ k 3 Tr(F 3 F 4 ) A 15 = k 3 -F x -F A -k 2 Tr(F 2 F 3 ) 
A14 = k 4 -Fx-F 3 -k 2 Tr(F 2 F 4 ) A X6 = k 4 ■ F 2 ■ F 3 ■ k 4 Tr(FiF 4 ) . 
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As an example application of the computer algorithms, one may check that the symmetri- 



sation operator of section 4.1 



5 4 b = (1 - 7112) (1 - 7r 34 )(l + 7ri 3 7r 2 4)(l + 7r 23 )(l - 7r 24 ) , 

cicts ELS 

5 4B ^i = 8A1 + ±A 2 - 4A 3 + 4A 4 + 8A 5 + 16A 6 

5 4B ^i6 = -6Ai + 6A 3 - 6A 5 - 12A 6 + |A 7 + 3A 8 + §A> + 3A 10 + §A U + 3Ai 2 
and has rank four. 
A. 2 Four fermions 

In dealing with the spinor wavefunctions Ui one has to face two issues: Fierz identities, and 
the Dirac equation. Fierz identities not only allow one to change the order of the spinors 
but also give rise to relations between different expressions in one spinor order. The Dirac 
equation often simplifies terms with momenta contracted into (ui'j^Uj) bilinears. 

In this section it is shown how to construct bases for terms of the form (k 2 or k 4 ) x 
U1U2U3U4. A significant simplification comes from noting that the Dirac equation allows 
one to rewrite (ui'y^Uj) bilinears into terms with lower n if more than one momentum is 
contracted into the 7^ . A good first step is therefore to disregard the momenta temporarily 
and find all independent scalars and two-index tensors built from m, ... ,114. From the 
SO(10) representation content, 

(5+)® 4 = 2- l + 6- B + 3- m + (tensors with rank > 2) , 

one expects two scalars and nine 2-tensors. The scalars are easily found by considering, as 



in 21 



Ti(1234) = (ui7 a n 2 )(u 3 7aw 4 ), 
T 3 (1234) = ( Ull abc u 2 )(u 3labcU ^ . 

and similarly for the other two inequivalent orders of the four spinors. (Note there is no T5 
because of self-duality of the 7®.) From Fierz transformations, one learns that all T3 terms 
can be reduced to T x by T 3 (1234) = -12T X (1234) - 24Ti(1324) and permutations, and the 
identity {Ha){ a p{^ a )~i)5 = implies that Ti(1234) + Ti(1324) + 7i(1423) = 0, leaving for 
example Xi(1234) and Ti(1324) as independent scalars. 

Generalising this approach to two-index tensors, it turns out that it is sufficient to 
start with 

Tu(1234) = (n l7 m n 2 )(n 3 7 n n 4 ), 
T 31 (1234) = (ui7 a 7 m 7 n ^)(u 3 7a«4), 
T 33 (1234) = (u 1 7 ab 7 m n 2 )(u 3 7a67^4), 
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and permutations of the spinor labels. It would be very tiresome to systematically apply 
a variety of Fierz transformations by hand and to find an independent set. Fortunately, 
by choosing a gamma matrix representation (such as the one listed in appendix |B|) and 
reducing all expressions to polynomials in the independent spinor components uj, . . . , uj 6 , 
this problem can be solved with computer help. As expected, one finds that the Tij(abcd) 
span a nine-dimensional space, and a basis can be chosen as 

Tn (1234), Tii (1324), T n (1423), T n (3412), T u (2413), T u (2314), 

T 31 (1234), T 31 (1324), T 31 (2314) . (A.l) 

A typical relation reducing the other Tij(abcd) to this basis is 

T 3 i(3412) = 2Tn(1234) - 2T n (3412) + T 31 (1324) + T 31 (2314) + 2^^(1234) . (A.2) 

Having solved the first step, it is now easy to include the two or four momenta, taking 
the Dirac equation into account. Consider first the case of two momenta. Starting from 
the two-tensors in (|A.1|), one gets the three independent scalars 



(«i|£ 3 ?/ 2 )(m 3 #iU 4 ) , (uih u 3)( u 2h u 4:) , (ui^ua,) (u 2 ^u 3 ) . 
In addition, there are four products of the two independent scalars T(1234) and T(1324) 



with the two independent momentum invariants si 2 and si 3 . By contracting ( |A.2j ) with 
momenta, one can show that 

si 2 Ti(1324) -si 3 Ti(1234) 

= -(■ui^ 3 ti 2 )(n 3 ^ 1 n 4 ) + {u l l/^ 2 uz){u2$ 1 u i ) - («i# 2 «4)(«2#i«3) , (A.3) 

and this relation can be used to eliminate si 2 Ti(1324). (It will become clear later that 

there are no independent other relations like this one.) There are thus six independent 
k 2 u\ ■ ■ ■ -U4 scalars: 

(ui^ 3 u 2 )(u 3 ^ 1 n 4 ) si 2 Ti(1234) 

(^ 2 u 3 )(u 2 ^u 4 ) si 3 Ti(1234) (A.4) 

(U1# 2 M4)(U2#1^3) si 3 Tl(1324) 

Note that there is only one completely antisymmetric combination of those, given by the 
right hand side of ( |A.3| ) . Similarly, in the case of four momenta, one finds ten independent 
k A u\ ■ ■ ■ U4 scalars: 

B 1 = 8X2 (lil#3 u 2)( u 3#l«4) B 2 = S i3 («l^ 3 U 2 ) {u^ x U A ) 

B 3 = S12 (U1#2'U3)(U2#1'«4) B 4 = Si 3 {uift 2 U 3 ) («2#1«4) 

B 5 = S12 (u^ui) (112^113) B 6 = si 3 (ni^ 2 n 4 )(u 2 ^ 1 ti 3 ) (A.5) 

S 7 = sf 2 Ti(1234) S 8 = si 2S i 3 Ti(1234) 

S 9 = s? 3 Ti(1234) Sio = s? 3 Ti(1324) 
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Working in a gamma matrix representation, it is again simple to construct a computer 
algorithm which reduces any given k 2 u± ■ ■ ■ U4 or k A u\ ■ ■ ■ 114 scalar into polynomials of the 
spinor and momentum components. The Dirac equation can then be solved by breaking 
up the sixteen-component spinors U{ into eight-dimensional chiral spinors u\ and u|, as in 
eq. QB.1| ), One obtains polynomials in the momentum components kf and the independent 
spinor components (u?) 1 '" 8 . However, a great disadvantage of this procedure is that it 
breaks manifest Lorentz invariance. For example, one encounters expressions which contain 
subsets of terms proportional to the square of a single momentum and are therefore equal 
to zero, but it is difficult to recognise this with a simple algorithm. The easiest solution 
is to choose several sets of particular vectors k{ satisfying kf = and ^ k{ = and to 
evaluate all expressions on these vectors. (By choosing integer arithmetic, one easily avoids 
issues of numerical accuracy.) Substituting these sets of momentum vectors in the bases 
( |A.4| ) and (A. 5) gives full rank six and ten respectively, showing they are indeed linearly 



independent. 

Equipped with a computer algorithm for these basis decompositions, one finds, for 
example, that the symmetriser 5 4 f of section 



5 4F = 4(1 + 7Ti 2 ) (1 - 7T 2 3)(1 + 7T24) , 

acts on the B\ . . . B\q basis as 

S AF B X = -125 4 + 125 5 + 12B 6 , 

S 4F B 10 = 8B 1 + 16B 2 - 8B 3 - 16B A + 8B 5 + 16B 6 - 2AB 7 - 24B 8 - 2AB 9 
and has rank three. 
A. 3 Two bosons, two fermions 

The combined methods of the last two sections can easily be extended to the mixed case 



of two bosons and two fermions. In the one-loop calculation of section 3.3, one encounters 
scalars of the form kui U2F2F4. A basis of such objects is given by 



(A.6) 



There are two combinations antisymmetric in [12] and symmetric in (34): 

-Ci + 4C 2 + C 4 and C 2 + C 3 . 

Finally, there are ten independent scalars of the form k 3 uiU2F^F4 (relevant to the two-loop 
calculation of section |4.3j ) , and they can all be obtained by multiplication of C\ . . . C5 with 
the two momentum invariants S12 and S13. 



Ci = 


(u l7 a u 2 )k 3 a FlF b 4 c 


c 2 = 


{ Ull a u 2 )Fl b Ftkl 


c 3 = 


{u^u^Flkl 


C4 = 


{ Ull abc u 2 )Fl b F' cd kl 


c 5 = 


( Ul r bc u 2 )F* b F* d k l d 
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B. A gamma matrix representation 

A convenient representation of the S0(l,9) gamma matrices is given by the 32 x 32 matrices 

(7 a ) Q/3 \ 



where 



■pa 

(/)«/? 
( 7 9)«/3 



(7°) 



ad 







ll6 = (7°)a/3 
(~ l80 ) = 

lJ 



-(7 9 



) a/3 



and (j a ) a P = — (^ a ) a p, a = 1 ... 8, is a real, symmetric 16 x 16 representation for the S0(8) 
Clifford algebra, 

a"' 



a/3 



(7' 



a = 1 . . . ; 



as given in appendix 5.B of [21]. The matrices r a satisfy the S0(l,9) Clifford algebra 
relations, 

{r a ,r b } = 2^132, v ab = (+ ), 

and bilinears of chiral spinors (with, say, positive chirality) are constructed as 

( nr [ai...a fc ] v) = ( n7 K...a fc ]^ = U <* ^) ^^fl . . . ( 7 a *% V S . 

This representation is particularly suitable for the calculations outlined in appendix [A| 
because it allows a simple decomposition of SO(l,9) spinors into SO(8) spinors due to its 
block structure: 

/is \ 



r° . . . r 9 



lie 

-i 16j 



r 1 • • • r s 



o -i 8 o o 
o o i 8 o 

\ -1 8 J 

Therefore, the Dirac equation for a chiral 16-component spinor u, 

can be solved by splitting u into two chiral eight-component spinors of SO (8), 



it 



u 



with 7 ' 



u 
u° 



—u c 



One obtains the coupled equations 

(d + d 9 )u s - (a ■ 3)u c = 
(d - d 9 )u c - (a T ■ d)u s = 
(with eight-dimensional dot products). These can be solved for u s in terms of u c 



(a ■ d)u c 



1 



-(a ■ k)u c , 



(B.l) 



where k + = —id + = ^7§(<% + dg). 
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